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Abstract 

We use a form of the fluctuation-dissipation theorem to derive formulas giv- 
ing the rate of production of spin- 1/2 baryons in terms of the fluctuations 
of either meson or quark fields. The most general formulas do not assume 
thermal or chemical equilibrium. When evaluated in a thermal ensemble we 
find equilibration times on the order of 10 fm/c near the critical temperature 
in QCD. 

PACS: 12.38.Mh, 25.75.Dw, 24.85. +p 



1 



Typeset using REVTgX 



I. INTRODUCTION 



Stable or long-lived baryons are readily detectable in high energy nuclear collisions, such 
as those at the CERN SPS (Super Proton Synchrotron) and at the Brookhaven National 
Laboratory RHIC (Relativistic Heavy Ion Collider). They have masses on the order of 1 
GeV, which makes them relatively sensitive to the temperatures and expansion rates in such 
collisions. Recent measurements with Au+Au collisions at y/s = 200A GeV at RHIC are 
consistent with protons, lambdas, sigmas, cascades, and omegas, and their anti-particles, all 
being in chemical equilibrium at a temperature of 170 ± 10 MeV jrj. This temperature is 
close to the expected critical or rapid crossover temperature in QCD between a quark-gluon 
phase and a hadronic phase. How is it possible to understand such near perfect equilibration 
on the relatively short time scales of high energy nuclear collisions? 

One approach is kinetic theory, which was reviewed in some detail already by Koch, 
Miiller and Rafelski ||. At finite temperatures, but still in the hadronic phase, there exist 
many mesons, especially tt, K, p and u. Two body reactions like tttt — > BB dominate at 
moderate temperatures, say 100 MeV and less. As the temperature goes up so does the 
density of mesons. Multi-particle reactions, such as rnnrKK — > BB become increasingly 
important. These multi-particle reactions cannot be measured directly in the laboratory, 
of course, but they can be related to the inverse reactions by detailed balance. Some of 
these inverse reactions, those involving proton plus anti-proton annihilation into mesons, 
were measured at LEAR (Low Energy Anti-proton Ring) at CERN. Rapp and Shuryak || 
have estimated that the sum of all reactions of the type nir — > pp are able to yield fast 
equilibration times for anti-protons at temperatures of order 150 to 170 MeV, perhaps as 
short as several fm/c. However, there are at least three difficulties with applying kinetic 
theory, using vacuum reactions rates, to high energy density matter. First, experimental 
data on the annihilation of hyperons into mesons is practically nonexistent. One must 
rely on approximate SU(3) flavor symmetry to estimate the rates based on nucleon data. 
Second, there is no sound practical method to implement microscopic reactions involving 
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initial states with more than two particles in a microscopic cascade computer code. Two 
body initial state reactions are assumed to occur when the colliding particles are within a 
distance \J cr/n of each other, where a is the cross section. For three or more particles the 
criterion becomes ambiguous. Third, and perhaps most importantly, the very concept of 
localized interactions occurring as in vacuum is no longer applicable. When particle densities 
reach 1 per cubic fermi with particles whose physical extent is also about 1 fm, it is not 
possible to define in and out states as in the vacuum. The interaction radius is comparable 
to the physical matter radius of the hadrons. The best that kinetic theory can do for hadrons 
at high temperatures is to tell whether the equilibration time is small or large but probably 
cannot give a quantitative number. 

Another approach is with DCC (disoriented chiral condensates) In this approach 
baryons are described as topological defects of the chiral field using the model of Skyrme. 
Domains are formed at some early time, and as the matter expands all the domains must 
eventually line up with the surrounding vacuum. During this process defects are formed with 
a probability that was first calculated in the context of the early universe || . The probability 
of defect production is inversely proportional to the cube of the domain size. Smaller domains 
yield more baryons and anti-baryons. Together with Wong, one of us showed that the 
anomalously large number of Q and fl observed at the SPS could be understood in terms 
of this mechanism |J. The typical domain size needed is 2 fm, which is just in the range 
predicted by many different approaches to DCC formation However, it is difficult to 
make much more quantitative calculations in this approach without very extensive numerical 
simulations. Even then, one may question whether a low energy effective meson field theory 
can reasonably describe variations over length scales as small as 2 fm. Plus the relevant 
frequencies involved are twice the proton mass. 

At first it would seem that these two approaches are almost orthogonal to each other. In 
fact they are trying to describe the same physics, baryon production at high energy density, 
starting from two opposite perspectives. Kinetic theory attempts to describe the process 
with many mesons that propagate freely between localized collisions. The DCC approach 
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assumes that the meson density is so high that they may be treated collectively as a classical 
field; baryons arise as topological defects of this field. 

In this paper we attempt a more general description than either of the above. We 
use a version of the fluctuation-dissipation theorem, basically the same formalism as used 
to compute the production rate of real and virtual photons in hot matter. The baryon 
production rate can be expressed in terms of the fluctuations in the local meson fields or 
in terms of fluctuations in the local quark fields. The formulas derived can be applied 
to systems in equilibrium or out of equilibrium. In the former case we evaluate the rates 
numerically; this is possible because of the Boltzmann weighting of states. In the latter 
case one must specify the probability of different states according to the preparation of the 
non-equilibrium system under consideration. We first compute the rates for non-strange 
baryons only, then we generalize to hyperons using SU(3) flavor symmetry in the interaction 
Lagrangian. 

II. NUCLEON PRODUCTION 

In this section we analyze the production of nucleons and anti-nucleons by vector and 
axial-vector fields or currents. We first derive general formulas that express the production 
rate in terms of fluctuations of these fields or currents. Formulas are given for arbitrary 
weighting of states. The formulas are then evaluated in thermal equilibrium with a Boltz- 
mann weighting of states. The vector and axial-vector fields or currents are first expressed 
explicitly in terms of pion fields using the nonlinear sigma model. They are alternately ex- 
pressed in terms of the quark fields using vector meson dominance to determine the relevant 
couplings. 

A. General formulas 

With a view toward the conventions used in the nonlinear sigma model we write the 
coupling of nucleons to vector V" and axial-vector A® currents or fields as 
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C = -^-^y« + g^rf^^Al , (1) 

where gr^ ~ 1.26 is the axial coupling constant relative to the vector. We do not need to 
specify anything more about V® and A a . They are given fields or currents that produce 
baryons via the above coupling. The production rate of nucleon/anti-nucleon pairs is calcu- 
lated in the same manner as dileptons M using linear response theory, equivalently a version 
of the fluctuation-dissipation theorem. We show explicitly the steps for production in the 
axial-vector channel. 

The matrix element for the process i — > /+ baryon + anti-baryon is 

S fi = g A (f\ J d 4 xA«(x)J 5 /(x)\i) . (2) 

Here Jf M is an abbreviation for the axial baryonic current as expressed in the Lagrangian 
above. The initial and final states are arbitrary, usually fully interacting states except for 
the specific interaction in Eq. (JJ). The axial baryonic current, corresponding to the baryons 
in the final state, is 



J5/{X) = T/V£^ e " lX ' (Pl+P2) ^(^'^)7V^b2, S2 ) • (3) 
As a consequence of the translational invariance we can write 

(/|A*(x)|i) = e fa -*(/|A-(0)|i>. (4) 

Here k = ki — kf is the difference of the four-momenta of the initial and final states. The 
transition rate is 
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Rfi = , (5) 

where T is the time interval and V is the volume. Explicitly 

~Y^UK{x + x')J^{x + x>)\i)(i\Al(x>)Jt(x>)\f) 
= ^Ihk I rf4xe "" (fc ~ Pl " P2) ^l^(°)K)(^'(°)l/) 

x w(pi, sihVy^fe, s 2 )v(p 2 , s 2 )l l v —u{p u Si) . (6) 



Averaging over initial states with an arbitrary weight Wi and summing over final states, we 
arrive at the following differential rate. 

dR = g^A E Wi J2(2nm Pl +P2 - k)(f\A«(0)\i) 

x (^(0)|/>Tr 
By introducing the correlation function 



h + m N 5 r a p 2 -m N 5 v r h 
2m N ' ' 2 2m N ' ' 2 , 



(7) 
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and evaluating the trace 
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we arrive at 



dR f 



N 



9\ 



h+m N 5 T a p 2 -m N 5 ^r 6 \ 

^ 7 7 M ir^ 7 7 17 

2m/v 2 2 ?7i 7\r 2 y 

P1V2 + P1P2 + - Pi • pa) 



(9) 
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■d 3 Pl d 3 p 2 2(2tt) 6 

Apart from the rather trivial Lorentz tensor coming from the form of the axial baryonic 
current, all the physics is contained in the correlation function. 

The corresponding rate arising from the vector interaction is easily obtained. The only 
differences are the replacement of the axial-vector with the axial field or current, A — > V, 
setting g\ = 1, and changing the sign of the mass-squared term in the Lorentz tensor. 
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dRv 
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V^ )aa ( P l + P2) [P^P U 2 + PiPt - <T (pi ■ P2 



m 



N 
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These are the first significant results in this section. They are quite general, but they require 
knowledge of the correlation functions for the initial states, appropriately weighted, and final 
states for the specific system and conditions under consideration. 

Now we evaluate the correlation functions in thermal equilibrium. The states \i) and |/) 
are conveniently assumed to be eigenstates of the Hamiltonian H with eigenvalues Ei and 
Ef. The weight is 
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where (5 is the inverse temperature and Z is the partition function. Following standard 
practice, the retarded correlation function in position space is 

A™*{x) = -i^Tr (e-^ [Al(x),Al(0)}) . (13) 

Its Fourier transform is 

AfJ a \p) = -i | o °° dx°d 3 xe^ p0+ ^ x °- x ^ 

x E {MiK(x)\f)(f\AUO)\i) - w f (i\Al(x)\f)(f\At(0)\i 

= -i J™ dx»d*xSv°- k0+i ^°-^-V\ J2w t (l- e^°) 



,/ 

x(<|A«(0)|/>(/|^(0)|<> 

w fl _ e Pk°\ ( 27r ) 3 ( 5 3 (p - k) 
= £ V po-fco + fe -<<|^(0)|/X/K(0)|<>. (14) 

As before, k = ki — kf. The imaginary part is 

W^p) = ~E^ (!" e/3p °) (^(p - fe)(<|^(0)|/>(/K(0)|<> . (15) 

Hence 
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Exactly the same relationship holds in the vector channel with the substitution of A with 
V. 

The sum of the rates of production by vector and axial-vector fields or currents can now 
be written as 

dR 1 1 

1 2 d 3 Pl d 3 p 2 ~ (27r) 6 e^+^) - 1 

x {im^! + pa) [p^ 2 + ftfi - gT (pi • p 2 + <)] 

+ 9 2 A lmA^ aa ( Pl + p 2 ) [rtfi + V \v» - gT ( P i • p 2 - rr&)] } . (17) 

This is the second significant result in this section. In order to make further progress, we 
need specific knowledge of the vector and axial- vector fields or currents and their correlation 
functions. 



B. Production in terms of pion fluctuations 



The self-interactions among pions in the nonlinear sigma model are contained in the 
Lagrangian 



where 



f 2 

C = ^Tr (dpUdPrf + mlU + m 2 ^ 



(18) 



u = e, 

e=exp Gr % 



(19) 
(20) 



Here ir a is the pion field and / T ~ 93 MeV is the pion decay constant. The pions couple to 
the baryons as in Eq. (^ with the derived currents 



K = -Tr 
" 2 



T d^ a + 

J TT 



ya 



-*Tr 
2 



_aoc_Qa _c i 

J TT 



(21) 
(22) 



To first order in the pion field, this gives the usual derivative coupling on account of the 
Goldberger-Treiman relation gx^N = f-wg-wNN- By expanding the interaction in powers of 
the pion field one may derive the contribution to the production rate involving two, three, 
four, five, etc. number of pions in the initial state. 



C. Production in terms of quark fluctuations 

Since pions are made of quarks it should be possible to express the vector and axial- 
vector fields or currents in terms of quark operators. We can make the connection using the 
hypothesis of vector meson dominance (VMD). In the two-flavor world the currents can be 
expressed in terms of the p, a\ and pion fields as 

V« = -g pN NPl, (23) 
= Wvw a? + g l! n Ld 

9A " 9A m N 
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Using the Goldberger-Trieman relation and VMD the quark currents are related to the 
hadronic fields. 

Q1, T Q = (25) 
T a vr? 

q^l 5 -q=^ L a a lll + pion (26) 

Apart from the pion pole we obtain the desired relationship. 

V; = "^f^y? (27) 

K = ^^Wy« (28) 

The imaginary part of the retarded correlator is given in terms of the spectral density p{s) 



as 



ImV^(k) = -3 [g, v - ^ 



n Pv (s = k 2 ), (29) 



with a similar expression for the axial channel. The factor of 3 arises from the sum over 
isospin indices. 

In our study of baryon/anti-baryon production, the spectral densities are needed only 
for s > 4:fn 2 N . Then s is large enough that they may be computed using perturbative QCD. 
(The vector one may be measured directly in electron-positron collisions; in the region of 
\/s of several GeV a sum of hadronic resonances gives essentially the same answer, a nice 
manifestation of duality.) They are U 



, x { gpNNgpTrnY s ( a s (s) \ 

pv{s)= {-^r) 8^1 1+ — ' (30) 

PA(s)=( 9 ^^) 2 ^(l + ^ + ..). (31) 

V JA< J 87T 2 V 7T J 

Here a s (s) is the QCD coupling evaluated at the scale s. 

According to KSFR |10[ 2g p f 2 = mi. Furthermore it is usually assumed that g p is 
universal in the sense that g P NN = gpnn, a result that holds rather well numerically. Then 
the coefficient of the vector spectral density is especially simple. 



In the absence of any better information it is quite reasonable to assume that the same holds 
in the axial-vector channel. 

l < ) "4/4 (33) 
This means that the vector and axial-vector currents have the same spectral density since 
the axial coupling g& then cancels out. Equal contributions is quite natural when the up 
and down quark masses are very small, as is the case in the real world. 

The tensor algebra can now be done using k = p\ + p 2 - It yields the third significant 
result of this section. 

dR _3 1 s(s-m? N ) ( as(s) \ ( , 

1 2 d 3 Pl d 3 p 2 8(27r) 7 e^+^)-l /4 [ l ^ n ^'") ^> 

This result is quite remarkable in that the rate is inversely proportional to the fourth power 
of fa and does not depend on any other hadronic parameters except the nucleon mass. 

Finally, the overall rate of production may be computed by integrating over the momenta 
of the outgoing nucleons. Since the threshold energy, 2m n, is much greater than the tem- 
peratures envisioned, T < 200 MeV, it suffices to drop the minus one in the Bose-Einstein 
distribution factor in the differential rate. Then, with K n denoting the Bessel function of 
the second kind, we get 

9 / a s (Am 2 N ) \ m%T 4 \m% 2 / m N 



R(NN) = — — 1 + ^ + 



(2tt)5 V J ft 



T 2 V T 

2 



(35) 



In the nonrelativistic limit this becomes 



This is the total rate for the production of pp, pn, np and nn. The individual rates are 
related as: R{pn) = R{np) = 2R{pp) = 2R(nn). 
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It should be noted that we have evaluated a s at threshold for nucleon/anti-nucleon 
production where the rate is a maximum. According to the latest analysis [IT] a s {m 2 T ) = 
0.35 ± 0.03. Since m T = 1777 MeV is very close to 2m^ we use that as the numerical value 
in our later calculations. The first perturbative correction to the spectral density is only 
about 10%. 

It should also be noted that we have included the interactions involving the isospin 1 = 1 
currents only. If the isospin 1 = currents were included too one might expect the rates for 
pp and fin to increase while the rates for pn and np to be unchanged. Indeed, this is what 
happens, and in fact the rates for pp and pn become approximately equal, but we defer the 
actual analysis to the next section. 

Finally, we point out that since nucleons are composite objects they have form factors. 
These form factors will multiply the above rates and will serve to decrease them to some 
degree. We will defer the determination of the form factors to the next chapter. The reason 
is that there is accurate data on pp annihilations, and to do a precise analysis we must first 
include the isospin 1 = current. 

III. NUCLEON AND HYPERON PRODUCTION 

An amazing fact in heavy ion collisions at RHIC is that hyperons are produced in rel- 
atively great abundance. Measurements indicate that they are, for all practical purposes, 
in chemical equilibrium at a temperature of 170 ± 10 MeV. Therefore it behooves for us to 
analyze hyperon production. In addition, we now will include coupling to the isospin 1 = 
vector and axial-vector currents or fields too. We side-stepped that contribution in the last 
section for clarity of presentation and to avoid making phenomenological estimates of the 
relevant couplings: assumption of SU(3) invariance of the interactions will help to resolve 
that issue. We will consider fluctuations in the SU(2) sector only as well as fluctuations in 
the full SU(3) sector. The results are rather different. 
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A. General SU(3) invariant couplings 



Flavor SU(3) is not nearly as good a symmetry as SU(2). However, experience over many 
decades of research has shown that it is usually a very good reproduction of experimental 
data to put all the flavor symmetry breaking in the mass terms but to insist that the 
interactions be flavor symmetric. That is what we shall do too. 

The vector meson nonet is conventionally written as follows |T2|,0|. 

\ 



V 



V 



P~ 

K*- 



P' 



K 



-2wg I uJs_ 



(37) 



The singlet and octet components are actually a mixture of the physical u and cf> mesons. 



u 8 = <t> cos dy + u sin dy 
u s — u cos 9y — <fi sin dy 



(38) 
(39) 



Ideal mixing occurs when the u has no ss component while the <fi is pure ss. The ideal 
mixing angle is tan#ideai = l/v2", or #id ea i ~ 35.3°. Experimentally the mixing angle seems 
to be about 39° [ TT|] . We will approximate the mixing as ideal to simplify formulas. Such 
fine details are not likely to be important in the context we have in mind, namely, heavy 
ion collisions. Therefore we use the nonet representation. 



V 



\ 



P~~ 

K*~ 



P_ , u_ 

V2 + V2 



K* + 
K*° 



\ 



K 



*o 



(40) 



A similar assumption of ideal mixing in the axial-vector meson nonet yields the following 
representation. 



A 



/ «£, A (1285) 

a i 



V2 + 



K° -A (1420) 



/i(1285) 



(41) 
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Finally we give the matrix of the baryon octet. 



B 



V2 ~*~ VE 



V 



-0 



P 

n 

-2A° 

/ 



(42) 



As is well known there are three types of SU(3) invariant couplings: the F and D types, 
so-called because they involve the correspondingly labeled group structure constants, and 
the singlet coupling, which involves the trace of the meson matrix. 



C 



mt 



= M» (1 - av ) Tl (B^ [V m , B}) + p v Tr (B>fB) Tr (V M ) 

+ a v Tr (BY {V m , B}) + g A a A TT (B7V {A„ B}) 

+ g A (l - a A )Tr (B^j 5 [A„ B)) + ^/3 A Tr (i^ 5 #) Tr (A, 



(43) 



The overall normalization of this interaction Lagrangian is determined by the coupling of 
the nucleons; see Eqs. ([!]) and (p3|). Four parameters are introduced: ay and a^, which 
determine the relative contributions of the D and F type couplings in the vector and axial- 
vector channels, respectively, and j3y and /3a, which determine the corresponding singlet 
contributions. 

There was evidence already in the 1960's that a a was about 2/3 |I2"|| . This has been 
confirmed repeatedly over the years. For example, in their analysis of the spin content of 
the nucleon Close and Roberts jnj] determined that a ~ 0.635. As another example, Klingl, 
Kaiser and Weise [[DJ use vector meson dominance together with SU(3) symmetry to deduce 
a. a = 0.68. We shall therefore fix = 2/3. 

The value of ay is determined by the requirement that the coupling of the <fi vector 
meson, which has already been taken to be a pure ss state, to nucleons vanishes: g^NN = 
T5fl . This is just one aspect of the OZI rule This requirement fixes ay = (1 — /3y)/2. 

We enforce the standard ratio of coupling constants of u and p vector mesons to nucleons, 
9uNN — 3g P NN, as follows from the quark model and the conventional definition of p and 
uj currents, see Eqs. (03) and fl45|). Thus, we determine (3y = 1. This further implies that 
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ay = 0. Finally, we require that the coupling of the nucleon to the /i(1420) meson vanishes, 
in analogy to the vanishing coupling of the <p meson to the nucleon. This condition fixes 
I3 A = (1 - 2a A ) = -1/3. 

The relative couplings in the vector channel with the choice ay = and (3y = 1 are shown 
in Table I. The absolute normalization may be inferred from the nucleon-nucleon couplings. 
The corresponding couplings in the axial- vector channel with the choice a a = 2/3 and 
(3a = — 1/3 are shown in Table II. 



B. Rates 

The invariant differential rates for all baryons in the octet can now be inferred. The 
only missing pieces are the spectral densities in the various channels. The currents as 
conventionally defined are, for example, 

$ = \(*fu-*fd) (44) 

3^ = \{uYu + d^d) (45) 
o 

J? = -\^s (46) 
and so on. The spectral densities for these currents from perturbative QCD are 

M.>-£(i + *a) W 

pjs) = >(l + «M) (48) 



72vr 2 V" ' 7T 
s / ' . at. (si 



(49) 



367T 2 V 7T I 

These are applicable above some threshold value which is always above the threshold for 
production of the corresponding baryon/anti-baryon pair, typically 1.5 to 2.5 GeV 2 . These 
threshold values could be estimated rather well by making use of the QCD sum rules ||. 

The rates are now determined from the SU(3) couplings, given in Tables I and II, and 
the spectral densities, given above. The relative weights are given in Tables III and IV from 
which the rates may be inferred. For convenience define the pair of functions 
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. 2 ^ann( s ) 2s 2 -(m 2 + m 2 )s-(m 2 -m 2 2 ) 2 ±Qm l m 2 s ( a s \ 

r±(mi ' m2) = ft I 1 + ~) 



(50) 



where the ± corresponds to vector /axial- vector contributions. The function Fann(s) is a 
form factor, alluded to in the previous section, and determined in the following subsection. 
Using the symbol r as shorthand notation for ExE 2 dRj d 3 p\d 3 p 2 some examples are given 
below. 

r(np) = 2r + (m N ,m N ) + 2r_(m N ,m N ) (51) 

82 

r(pp) = 2r+(m N ,m N ) + —r-(m N ,m N ) (52) 

ol 

25 

r(Ap) = 3r + (rn A ,m N ) + —r-(rn A ,m N ) (53) 
r(H~A) = 3r + (mj v ,m s ) H — -r_(mA,m=) (54) 

Altogether there are 46 nonvanishing combinations of baryon/anti-baryon pairs. 

The rates cannot in general be evaluated in closed form with the form factor included. 
However, it turns out that a very good approximation (within 10% at T = 200 MeV) is to 
evaluate -^ann( s ) a ^ the average value s = (mi + m 2 ) 2 + 3(mi + m 2 )T, as discussed in the 
next subsection. Let us define 

d 3 pi d 3 p 2 



R 



± 



/ -p F~ T± ■ (55) 

J rj\ J2/2 



If we can evaluate F| NN at the average value of s then the integral can be done in closed 
form. It is 

R± = 9(1 j7^f4 )r8 ^2 {^K^K^z,) + iziK^K^) 

4 ( 27r ) Jn 

±( Zl ± z 2 ) 2 K 1 (z 1 )K 1 (z 2 ) + [16 + ( Zl ± z 2 ) 2 }K 2 ( Zl )K 2 (z 2 )} F 2 NN (s) , (56) 
where Z{ = rrii/T. 
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C. Form factors 



Nucleons are composite objects, hence they have form factors that depend on the specific 
process. The appropriate form factor here is not the electric or magnetic form factor since 
we are coupling the nucleons to mesonic currents, ultimately expressed in terms of the quark 
fields. We can obtain a very good estimate of the relevant form factor by comparing the 
rate for pp production, as derived above, with the rate for pp annihilation, as obtained from 
kinetic theory. 

The kinetic theory expression for the annihilation rate is 



E ' E2 d^¥p- 2 {PP) = (2^ /(El)/ ™ <NN(S) ' (57) 

where cr A ^ N (s) is the annihilation cross section which specifically excludes a baryon/anti- 
baryon pair in the final state. If we approximate the thermal distributions f(E) = 
exp(—E/T), which is a very good approximation to the Fermi-Dirac distribution at the 
modest temperatures of relevance here, we obtain the following simple expression. 

g^-Rann /_ n 2 



ElE2 d*p7w PT>) = (2^ exp ^ El + E *)l T ^< s - H)^ann(^) (58) 
In chemical equilibrium the rate for production must be equal to the rate for annihilation. 
The former is given in Eq. (|52|) . Dropping the f\ contribution for simplicity of presentation, 
meaning that the factor 82/81 is set to 1, which is an approximation to the rate better than 
0.5%, we get 

EA l^k^ - 32(W eXP '-< £ ' + E ^Ji i 1 + T) F -" W ' (59) 
Equating these two yields an expression for the form factor in terms of the annihilation cross 
section. 



The experimentally measured values of the pp annihilation cross section from several hundred 
MeV/c to 8 GeV/c lab momentum has been meticulously parameterized by Cugnon and 
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Vandermeulen see also the review by Dover, Gutsche, Maruyama and Faessler [18 . 

The fit is 

^ann(Pl) = -7= + ;tt (61) 



given in mb when the lab momentum pi is given in GeV/c. 
We have fit the form factor with the monopole function 

Fann(s) = 2.21 + (s- Am%) /A* ' (62) 

where A = 1.63 GeV. This function gives a very good representation for s — Am 2 N > 0.5 
GeV 2 , but overestimates Fann by about 10% at s — 4m 2 N . This overestimate is acceptable 
because, as we shall see below, the average value of s — <im 2 N is greater than 0.5 GeV 2 for 
T > 100 MeV. 

For hyperons we choose the parameterization of the form factor to be 

Fann(s) = 2.21 + [ S -(m 1+ m 2 ) 2 ]/A 2 ' (63) 

with the same value of A. 

The total and elastic cross sections for pn have been measured, and the annihilation cross 
section estimated, for lab kinetic energies between 450 and 1068 MeV [19[] and at 3.5 GeV 



BOfl. The annihilation cross section for np has been explicitly measured for lab momenta 



between 100 and 500 MeV/c [21]. In all these cases the annihilation cross section for pn 
and np has been equal to the annihilation cross section for pp, albeit with large error bars 
in the two former cases. The difference in the thermal production rates between pp and np 
calculated here is well within the error bars. [Note that the rates in Eqs. (|5T| ) and ([52] ) do 
not differ much because r_(m7v, v^n) is rather small compared to r + (mjy, m N ) for s > ^m 2 N ] 
This is a gratifying conclusion which points to the consistency of our results compared to 
experimental data. 

In principle the differential rates must be integrated over all energies, including the s 
dependence of the form factor. A rough approximation is to evaluate the form factor at the 
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thermal average value of s, which is s = (mi + m-i) 2 + 3(mi + m-ifT, in the nonrelativistic 
limit and dropping terms of relative order (T/m) 2 . It turns out that this approximation is 
good to better than 10% for temperatures less than 200 MeV. 



IV. NUMERICAL RESULTS 

Consider a system at fixed temperature and volume but not necessarily in chemical 
equilibrium with respect to the baryons. If one of the baryons, say the anti-proton for 
definiteness, is out of chemical equilibrium for any reason, how long does it take for it to 
return to equilibrium? The rate equation for the density is 



_ n p n b 



equil equil 



(64) 



where np qml is the equilibrium density for anti-protons and n^ qml is the equilibrium density 
for the baryon species b. The characteristic time scale for bringing the anti-protons to 
equilibrium is 

r p = n^/J2R(bp). (65) 

b 

This characteristic time is more intuitive than the rates themselves. A fully dynamical model 
of the evolution of matter is required for the detailed knowledge of how the abundances de- 
velop, but for the purpose of gaining insight to the dynamics the characteristic equilibration 
time is perhaps more useful. 

We plot this time for the proton, lambda, sigma, and cascade baryons (same as for 
the anti-baryons in net baryon-free matter) in Figs. [I] and 0. Figure [l] shows the times 
when only fluctuations in the SU(2) meson sector (no strangeness) are allowed while Fig. |2| 
shows the times when all mesons or currents are included (including those with strangeness) . 
The equilibration times are strongly decreasing functions of increasing temperature. That 
is typical of thermal processes; thermal rates are generally strongly increasing functions 
of temperature because of the dominant Boltzmann factor. The time for nucleons is the 
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shortest in the SU(2) case, which is quite natural since nucleons contain no strange quarks 
and they are the lightest baryon species. The lambda has the longest equilibration time as 
a consequence of the magnitude of its couplings to the SU(2) fluctuations. Going to the full 
SU(3) fluctuations, see Fig. ^|, shortens the cascade time considerably but, surprisingly, the 
lambda now has the shortest equilibration time! It is interesting to note that the £77(3) 
symmetry is broken not only by the different masses of the baryons but also by the mixing 
of the singlet and octet mesons. 

As mentioned already, to compare with data from heavy ion collisions requires solving 
rate equations in an expanding and cooling system. Examples of how this may be done 
is described in P, p2|j23[| . In addition, feed down from the decay of higher mass baryon 
resonances will contribute to the observed yields. But to get a rough idea, suppose that 
the expansion time scale is about 10 fm/c and that there are full SU(3) fluctuations in the 
system. Draw a horizontal line at 10 fm/c in Fig. 0. The intersection with the various 
baryon species would suggest that these baryons would reflect a freezeout temperature in 
the range 168 to 180 MeV, the exact value depending on the species. This is approximately 
the range of chemical equilibration temperatures recently seen in Au-Au collisions at 130 
and 200 GeV at RHIC |TJ. 

V. CONCLUSION 

In this paper we have calculated the production of spin- 1/2 baryon/anti-baryon pairs 
through fluctuations in the strong interaction currents. The most basic formulation used 
a version of the fluctuation-dissipation theorem that does not rely on the system being in 
thermal equilibrium. If one has a model for these fluctuations those formulas may be used 
directly. We evaluated them in thermal equilibrium, which gives rise to equilibration times 
short enough that nucleons and hyperons may very well be in chemical equilibrium in heavy 
ion collisions at RHIC energies. 

Two natural extensions of our work arise. The first is to carry out the analogous cal- 
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culation for the spin-3/2 baryon decuplet. Coupling of the strong interaction currents to 
spin-3/2 baryons is much more uncertain than the coupling to spin- 1/2. The second is to 
apply the formulas derived in this paper to a dynamical model of the expanding matter. 
Only then will we be able to make direct contact with RHIC experiments. 
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Table I: Relative strength of vector couplings for ay = and (3y = 1- The names of vector 
mesons that couple to the given baryon/anti-baryon pairs are shown explicitly. 
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Table II: Relative strength of axial-vector couplings for ola = 2/3 and Pa = —1/3. The 
names of axial- vector mesons that couple to the given baryon/anti-baryon pairs are shown 
explicitly. The notation f\ stands for /i(1420). 
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Table III: Numerical vector channel multipliers in the expressions for the rates (ay = and 
Pv = 1). 
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Table IV: Numerical axial-vector channel multipliers in the expressions for the rates (a a = 
2/3 and ^ = -1/3). 
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FIGURES 
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FIG. 1. Numerical results for equilibration times in the case when only the fluctuations in the 
SU(2) meson sector are taken into accout. 
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FIG. 2. Numerical results for equilibration times in the case when all fluctuations in the SU(3) 
meson sector are taken into accout. 
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